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We study the flux and the angular power spectrum of gamma-rays produced by Dark Matter 
annihilations in the Milky Way (MW) and in extra-galactic halos. The annihilation signal receives 
contributions from: a) the smooth MW halo, b) resolved and unresolved substructures in the MW, 
c) external DM halos at all redshifts, including d) their substructures. Adopting a self-consistent 
description of local and extra-galactic substructures, we show that the annihilation flux from sub- 
structures in the MW dominates over all the other components for angles larger than 0(1) degrees 
from the Galactic Center, unless an extreme prescription is adopted for the substructures concen- 
tration. We also compute the angular power spectrum of gamma-ray anisotropies and find that, for 
an optimistic choice of the particle physics parameters, an interesting signature of DM annihilations 
could soon be discovered by the Fermi LAT satellite at low multipoles, I < 100, where the dominant 
contribution comes from MW substructures with mass M > 10^ A/©. For the substructures models 
we have adopted, we find that the contribution of extra-galactic annihilations is instead negligible 
at all scales. 

PACS numbers: 



I. INTRODUCTION 

Despite the compelling evidence for Dark Matter 
(DM), we know very little about its nature. It is com- 
monly assumed that DM is composed of Weakly Interact- 
ing Massive Particles (WIMPs) , that are kept in thermal 
and kinetic equilibrium with baryons in the early uni- 
verse through their weak coupling with ordinary matter, 
and that subsequently decouple from them when the self- 
annihilation rate drops below the expansion rate of the 
universe, thus naturally achieving the appropriate relic 
density. Being proportional to the square of the DM 
number density, the self-annihilation rate is today very 
small, but it can still lead to detectable signals in regions 
with very high DM density 

Indirect DM searches are based on the detection of par- 
ticles originating from DM annihilation or decay. If one 
focuses on the most widely discussed candidates, i.e. the 
supersymmetric ncutralino and the so-called B^^^ in the- 
ories with Universal Extra Dimensions, the mass of the 
DM particles should be approximately in the 100 GeV 
- 1 TeV range, and the characteristic energy of the pro- 
duced particles, roughly an order of magnitude smaller. 
At these energies, photons are among the best mcssen- 
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gers, since gamma-rays travel in the local universe along 
geodesies without significant energy losses. 

The most obvious target for indirect DM searches is the 
center of our Galaxy, where the large concentration of or- 
dinary matter is believed to be associated to a large DM 
density enhancement. However, the search for gamma- 
rays from DM annihilations in the Galactic Center (GC) 
is complicated by the presence of a strong signal of astro- 
physical origin, due to the point sources detected by the 
EGRET satellite d i, Q and by the H.E.S.S. telescope 
0. A further degree of complication is constituted by 
theoretical uncertainties on the DM profile in the inner- 
most regions of the MW halo which is inevitably affected 
by the presence of a SuperMassive Black Hole and by the 
surrounding distribution of stars 0, Q . 

To avoid such difficulties, different authors have sug- 
gested alternative strategies to unambiguously detect the 
signature of DM interpretation. One possibility is to look 
for enhancements in the gamma-ray flux with no obvious 
astrophysical origin. Moreover, the presence of possible 
distinctive features in the gamma-ray energy spectrum 
(as lines [1, [l3l or "bumps" [lH) may help in ruling out 
alternative, more conventional, interpretations. Ando et 
al. p^ . [isj have recently pointed out that unique an- 
nihilation features may also be detected in the angular 
correlation properties of the Extragalactic Gamma-ray 
Background (EGB). From an observational point of view 
the EGB is obtained by subtracting the galactic gamma- 
ray flux produced by cosmic-rays interacting with the 
interstellar medium at high galactic latitude from the to- 
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tal background measured by EGRET . The resulting 
EGB is isotropic and its energy spectrum is rather un- 
certain, especially at energies above 10 GcV [l5l[T6l. [l7l| . 
Unresolved extra-galactic gamma-ray sources like blazars 
certainly contribute to the EGB. However, current un- 
certainties on the luminosity function of these objects do 
not allow to exclude additional contributions from more 
exotic processes like DM annihilation within and outside 
our Galaxy. 

The annihilation signal along a given direction in the 
sky is contributed by four different sources: 

• the smooth MW halo, that is expected to con- 
tribute most towards the GC where the DM density 
enhancement is expected to be very large (see ref- 
erences above), 

• resolved and unresolved substructures within the 
MW halo (see [l^l and references therein). 

• extra-galatic DM halos at all redshifts [l^, includ- 
ing 

• their substructures. 

Ando et al. have shown that if DM annihilations in 
extra-galactic halos and their substructures contribute 
more than 30% to the EGB, the recently launched Fermi 
LAT satellite should be able to identify this 'exotic' con- 
tribution in the angular power spectrum. 

Their analysis does not account for the possible con- 
tribution of galactic substructures to the angular corre- 
lation signal [l2l . [l^ . Such contribution has been studied 
recently by Siegal-Gaskins who concluded that detectable 
features in the EGB angular power spectrum may exist 
and they can be used to infer the presence of galactic 
substuctures and constrain their abundance [2]| (in this 
case, annihilations in extra-galactic halos have been ne- 
glected). 

Our model for galactic substructures is taken from Ref. 
, which is inspired by the high resolution N-body sim- 
ulations aimed, precisely, at resolving subclumps within 
galaxy-sized halos (see Refs. jH,!^ ^^'^ the more recent 
Refs. [13, [lll,!!^])- At z = 0, substructures are predicted 
to span 12 orders of magnitude in mass. Thus, the de- 
scription of smaller objects is based on the extrapolation 
from the results at higher masses, since resolving such 
small clumps would be well beyond the resolution limit 
of the N-body experiments. In the attempt of limiting 
the theoretical uncertainties involved in this extrapola- 
tion, we have adopted two different prescriptions for the 
concentration of DM halos and subhalos of mass M as 
a function of redshift z, c{M,z). Then, to compute the 
gamma-ray fluxes from the resulting ~ 10^^ galactic sub- 
structures, we have used the hybrid analytical and Monte 
Carlo approach described in Ref. [T8| . 

The extra-galactic flux has been modeled using the for- 
malism of Ref. IT^ . adapted to our calculation setup. In 
this way, we were able to compute for the first time in 



a self-consistent way the normalization of the DM extra- 
galactic signal relative to the DM galactic foreground, a 
quantity that does not depend explicitly on the particle 
physics parameters (a weak implicit dependence remains, 
since the lower limit of integration of the subhalo mass 
function is sensitive to the particular particle physics sce- 
nario) . 

We have then computed the angular power spectra by 
running the public code HEALPix on our mock maps of 
galactic substructures, and comparing our result to those 
of Ref. OJ] and Ref. [2l|. 

The paper is organized as follows: in Sections HHIIVI 
we describe the theoretical setup, the properties of our 
model substructures and compute the mock full-sky maps 
of the gamma-ray flux due to DM annihilations, both 
galactic and extra-galactic. In Section|V|we compute the 
angular power spectrum of the mock gamma-ray sky and 
simulate its measurement using Fermi LAT in Section IVTl 
Finally, in Section |VlT] we draw our main conclusions and 
discuss our results. 



II. SIMULATION OF GALACTIC 
SUBSTRUCTURES 



The Cold Dark Matter scenario predicts the formation 
of a large number of DM halos, virialized structures with 
masses possibly as small as ~ 10~^ M© [S^, [2^. Such 
halos merge into larger and larger systems, leading to 
the formation of today's halos of galaxies and clusters 
of galaxies. A fraction of the small halos survives to 
dynamical interactions with the stellar and dark compo- 
nents until the present epoch. Earth-size substructures 
have indeed been found in numerical simulations within 
DM halos of mass 0.1 M© at a redshift of 74 [H. 

Resolving such small structures within a galaxy-sized 
halo today is out of computational reach, so that their 
spatial distribution, mass function and internal structure 
can only be estimated upon a rather uncertain extrapo- 
lation from the properties of higher mass satellites. Here, 
we work under the following hypotheses: 

• we assume that substructures trace the mass of the 
MW, i.e that their radial distribution follows the 
mass density profile of the parent halo (as found 



g. in Ref. [2J|, see however the discussion below 



on the results of the most recent numerical simula- 
tions), 

• wc assume that the substructures mass func- 
tion is well approximated by a power-law 
dn(M)/dln(M) cx normalized so that 10% of 

the MW mass lies in objects within the mass range 
IO-^A'/mw - IO-^Mmw, as found in Ref. [H. 



Under these assumptions, the number density of sub- 
halos per unit mass at a distance R from the GC can be 
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written as: 



Psh{M,R)=AM- 



0{R- 



(i?/rMW)(l + i?/rMW)2 



M0-lkpc-^ 



(1) 

where is the scale radius of our Galaxy and the cfTect 
of tidal disruption is accounted for by the Hcavisidc step 
function 9{R — rmin{M)). The tidal radius, rmin{M), 
is found by using the Roche criterion to decide whether 
a subhalo of mass M survives tidal interaction with the 
host halo (see Ref. According to our normalization 

about 53% of the MW mass is condensed within ~ 1.5 x 
10^^ subhalos with masses in the range 10~^ — 10^" M0. 
Their abundance in the solar neighborhood turns out to 
be ^ 100 subhalos pc""^. 

As far as the smooth component of the Galactic halo 
is concerned, we assume that the MW halo follows a 
Navarro, Frenk & White (NFW) profile [H, with 



a scale radius 



21.7 kpc and that the total mass 



enclosed in a radius r2oo, corresponding to the radius 
where the halo density is 200 times the critical density 
of the universe, is Mmw = 10^^ M©. An NFW fit to the 
DM halo of our Galaxy is consistent within 10% with the 
results of the Via Lactea I simulation [28l |. A popular 
way to characterize the mass distribution within the DM 
halo is using the so-called concentration, a shape param- 
eter defined as the ratio between the virial radius and the 
scale radius, c = r^Qo/fs- 

Subhalos are also assumed to follow an NFW profile, 
with a concentration that depends on their mass [2^ . 
To determine the dependence of concentration on the 
mass, we follow the prescription proposed in Ref. p9l | 
(hereafter BOl) in which the concentration parameter is 
found to depend on both the subhalo mass Mh and on 
its collapse redshift, Zc, defined as the epoch in which 
a mass scale Mh enters the non-linear regime. The col- 
lapse occurs when a{Mh)D{zc) ~ 1, where a{Mh) is the 
present linear theory amplitude of mass fluctuations on 
the scale Mh and D{zc) is the linear theory growth factor 
at the redshift Zc- In an attempt to bracket our theoret- 
ical uncertainties in modeling c{Mh,Zc) and its extrap- 
olation at low masses we have implemented two rather 
extreme models chosen among those considered in Ref. 
[l8l |. namely: 

• Bzo,ref- that extrapolates the c{Mh, Zc) relation of 
BOl below lO^'Mo with a simple power-law. 

• Bzf^ref'- which assumes that surviving subhalos 
do not change their density profile since their for- 
mation. Thus the concentration parameter at 
can be obtained from the one at z = through 
c{Mh,Zc) = c{Mh,z = 0)/(l + Zc). The values of 
c{Mh, z = 0) corresponds to those of the B^g^ref 
case and the collapse redshift Zc is obtained by ex- 
trapolating the expression proposed by BOl below 

WMq. 

Finally, the concentration parameters are not uniquely 
defined by the halo mass. Rather, they follow a log- 



normal distribution with dispersion Uc = 0.24 [29j and 
mean c(M): 



P(c(Af),, 



1 _/ln(c)-ln(e(JU» V 

e V ^^<'= / . (2) 
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A. Recent highlights from numerical simulations 

Two new sets of very high resolution numerical experi- 
ments have been recently released, namely the Via Lactea 
II [lOl and the Aquarius [3l|, simulations. The main 
characteristics of the MW subhalo population in these 
two simulations are summarized in Table H] and compared 
to those of our subhalos. Among the main differences, we 
note that the mass fraction in substructures is a factor 2-6 
smaller than that used in our model. Another difference 
is in the subhalo distribution which, in our case, trace 
the smooth mass distribution of the MW halo, pmwWj 
while Refs. [sH . [s^ suggest to use an Einasto profile (33j 
and in Via Lactea lithe innermost regions are best fitted 
by a power-law (1 + r)^-^ . 

The concentration parameter of subhalos in the two 
simulations depends from the distance from the GC and, 
in the Aquarius case, follows the prescription in Ref. 
(NOT) rather than the BOl model. Finally, in Aquar- 
ius, the subhalo density profile is also parametrized as 
an Einasto. rather than an NFW, profile. 

These differences, in particular the fact that neither 
in Via Lactea II nor in Aquarius the spatial distribution 
of subhalos trace the mass of the parent halo, do have 
a significant impact on the angular power spectrum of 
the gamma-ray flux, as we will show in an upcoming 
publication. Here, we perform all calculations under the 
assumptions listed above, and discuss how the differences 
among the subhalo models in Table U are expected to 
affect the angular power spectrum (Section I VII[) . 



III. MODELING THE GAMMA-RAY FLUX 

A. Gamma-ray flux from galactic subhalos 

The gamma-ray flux expected from the annihilation of 
DM particles can be written as: 

^{E„l,b) = ^(E,) X (/,6) (3) 



where the term 



X / 7 



(4) 



contains the dependence on particle physics parameters, 
while 

$'=°™°(V',/,6) = / dM [ dc [ [ ded<j) [ d\ 

JM Jc J JAn Jl.o.s 
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this paper 


Aquarius 


Via Lactea II 


dumpiness (< r2oo) 


53% 


8% 


26% 


Nsubhalos{< '"200 ) 


1.5 X lO^** 


2.3 X 10" 


7 X 10^^ 


dn/dM 




cx M-'-'> 




nsh{R) 


NFW 


Einasto 

QEinasto = 0.68 


cx {l + R)~^ 


c{M, z) 


BOl 


N07 


BOl 


pMwir) 


NFW 


Einasto 

QEinasto =0.21 


NFW 


PhaloiR) 


NFW 


Einasto 

QEinasto = 0.16 


NFW 



TABLE I: The main characteristics of the subhalo model con- 
sidered in this work compared with those measured in the 
recent Aquarius and Via Lactea II N-body simulations. The 
dumpiness is defined as the fraction of the dark mass in sub- 
structures within the virial radius of the MW r2oo = 210 kpc. 
Nsubhaios is the total number of substructures within r2oo. 
dn/dM is the subhalo mass function, nsh{r) their radial dis- 
tribution and c(M, z) is the model used for the concentration. 
Pmw and phaio represent the mass density profile of the host 
halo and its substructures, respectively. 



[pshiM,RiRQ,X,l,b,e,(t,)) X P(c)x 



X ^Zr:°iM,c,riX,l,b,0,4>)) X J{x,y,z\X,0,4>)] (5) 

represents the contribution to the foreground emission 
from the subhalo population and 

$™(M,c,r)= / / dcf,'de' f dX' 
J JAn Ji.o.s 



'pL,,{ M,c,riX,X'J,b,e'q^')) 
A2 



J{x,y,z\X',9'cl)') 



(6) 



is the contribution from a single subhalo. 

The total galactic flux is obtained by adding to Eq. [5] 
the contribution of DM annihilations in the smooth NFW 
halo of the MW. 

In Eq. m is the DM particle mass, CTanni' the an- 
nihilation cross section, and dNj^/dEj the differential 
photon spectrum per annihilation relative to the final 
state /, with branching ratio Bf that we take from Ref. 
[35| . In this work we adopt a rather optimistic particle 
physics scenario (in the sense that it provides large an- 
nihilation fluxes), in which the DM particle has a mass 
~ 40 GeV, a cross section cru = 3 x lO^^^cm'^s^^ and 
particles annihilate entirely in bb. 

In Eq. [5l Ail = 9.57 x lO^'^sr is the solid angle con- 
sidered for the integration, corresponding to the angular 
resolution of the Fermi LAT satellite, {l,b) are the galactic 
coordinates of the direction of observation, and J is the 
Jacobian determinant of the transformation between po- 
lar and cartesian coordinate systems. The galactocentric 



distance, R, can be written as a function of the coordi- 
nates inside the observation cone (A, 9, (f) and of b) 

through the relation R ~ \J + -^©^ ^ 2XRqC, where 
i?0 is the distance of the earth from the GC, and C is the 
cosinus of the angle between the direction of observation 
and the direction of the GC. 



B. Gamma-ray flux from extra-galactic structures 

Aside from local DM, the gamma-ray flux in any given 
direction receives a contribution from all the structures 
at all redshift along the line of sight. Adapting the for- 
malism of Ref. [l^ , we estimate here the contribution of 
extra-galactic structures, including the presence of sub- 



structures following the B 



and B 



models. 



Following Ref. [19|, the infinitesimal volume dV at a 
redshift z can be written as 



dV 



R^r^drdn 

(1 + Z)3 ' 



(7) 



where dQ is the solid angle, dr is the infinitesimal comov- 
ing depth and Rq is the scale factor at the present epoch. 
If, for a moment, we neglect the presence of subhalos 
and assume that the gamma-ray emission is isotropic, we 
can compute the number dN-y of gamma-ray photons pro- 
duced in dV in a time interval dt with an energy between 
E and E + dE and collected by a detector with effective 
area dA by integrating the single halo emissivity over the 
halo mass function j^{M, z): 



dN^ 



-t{z.E„) 



(l + z)= 



dM-^{M,z) (8) 



dAf. 



dVdA 



where _Bo and dt^ are, respectively, the energy and the 
time interval over which the photons are detected on 
earth. These quantities are related to those at the red- 
shift of emission through Eo ~ E/{1 z) and dt^ = 
(1 4- z)dt, so that dt^dE^ — dtdE. The halo mass func- 
tion, dn/dM, represents the comoving number density of 
DM halos of mass M at redshift z and the factor (1 -I- z)^ 
converts comoving into physical volumes. In the Press- 
Schechter formalism [301. the halo mass function is 



dlogM' 



(9) 



where p^, is the critical density, SIo.jti is the mass density 
parameter, v = ^^[jy-, cr{M) is the rms density fluctua- 
tion on the mass scale M and 5sc represents the critical 
density for spherical collapse. We refer to Refs. [l^, [s^l 
and for the exact computation of these quantities. 

dAf-y/dE represents the number of photons with en- 
ergy between E and E + dE produced in a halo of mass 
M at redshift z. The exponential e"'^'^^'^") is an ab- 
sorption coefficient that accounts for pair production due 
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to the interaction of the gamma-ray photons with the 
extra-galactic background Hght in the optical and in- 
frared bands. Following Ref. we adopt here the 
following expression that accounts for current observa- 
tional constraints: t{z,Eo) z/[3.3(£'o/10 GeV)~°-^] . 
The mean flux is obtained by integrating Eq. [8] along the 
line-of-sight: 



(10) 



dEodn I ' ' dEodtodAdn 

/ drRoe-^^''^°^ 
47r J 

jdM^{M,z)^{E,{l + zlM,z) 

c f , e-^(^'-^«) 

— / dz -— 

Att J Hoh{z) 

dM-^{M,z)^{Eo{l + z),M,z), 

where the last expression has been obtained by trans- 
forming comoving distances r into redshifts z, through 
the introduction of the Hubble parameter Hoh{z) = 
Ho\/^o,m{i + z)-^ + rio.A, where Hq is the Hubble con- 
stant and r^o,i the abundance in units of the critical den- 
sity at z = 0. 

The number of photons emitted in a single halo, 
dAf.y/dE, depends on the DM density profile (NFW in 
our case) and on the particle physics scenario (particle 
mass m^, annihilation cross section av and differential 
energy spectrum per annihilation. dN^/dEo). The NFW 
profile of a halo with mass M is completely specified by 
the concentration parameter and the virial overdensity 
l^vir of the halo. The virial radius r^ij., as opposed to 
the previously defined r200 7 is the radius of the sphere 
which enclosed an average density A„ir x Pm.. To be con- 
sistent with the notation used in Ref. [131 we define the 
concentration parameter as c(M, z) = rmr/i's and change 
the corresponding values in models Bz^^ref and Bz^.ref 
accordingly. 

Bearing all this in mind, we can express dM-y/dE as 

dM^ av dN^{E) M AyirPcr^miz) 

-{E,M,z) = — — (11) 



dE 



2 E ml 3 



c^{M,z) 



T ( 11^^ r2/2(a;„,.n,c(M,z)). 

/i(a;„,.„,c(M, z)y 
In the previous expression, the virial overdensity is [l^ : 
IStt^ + 82(f7„,(z) - 1) - 39(f7,„(z) - if 



(12) 

and the integrals Ii and I2 have an analytic expression: 



/„(a;„i„,a;,„,,) 



where g{x) — x "'^(1 -f- a;) 



g^-x dx, 



(13) 



In Eq. Ill[ the lower integration limit is set at the min- 
imum radius within which the annihilation rate equals 
the dynamical time: x^i^ ~ 10~^kpc/rs (?'« is the scale 
radius in kpc). We have checked that the results are not 
sensitive to a different choice for a;„,i„. 

Putting Eqs. [TTjandfTOltogether. we obtain the expres- 
sion for the isotropic gamma-ray fiux from extra-galactic 
DM halos: 



d<P 



dEodn 



(Eo) 



av c Pcr^o. 



SttHo 



(14) 



dz{l + z) 



, A^{z)dNy{Eo{l + z)) 
h{z) dE 



with 



A^(z) = / dM 



,v{z, M)f{v{z,M)) 
a{M) 



da 



dM 



AlAz,M) 



(15) 



and 



Al,(z,M) = J dc'P(c(M,z),c')^||^(c')W. 

] (16) 

A^^(z, M) represents the enhancement in the gamma-ray 
fiux due to the presence of a DM halo with a mass M at a 
redshift z. In A^(z) all these contributions are integrated 
over the halo mass function. Therefore A^(z) quantifies 
how much the annihilation signal is boosted up by the 
existence of virialized DM halos. 

To account for the presence of substructures in extra- 
galactic halos, we assume that a given fraction of the 
halo mass is concentrated in substructures with the same 
properties as the galactic subhalos described in Section 
im For consistency with the MW case, we require 10% 
of the parent halo mass M to be in subhalos in the mass 
range 10~^M — 10~^M. Following this requirement the 
mass fraction f{M) in substructures within a host halo 
of mass M can be fitted as 



/(A/) 



M 



logio(A//MQ) 
30 



(17) 



With the above definition, the presence of substruc- 
tures is taken into account replacing A\j in Eq. 1161 with 
the following expression: 



Al,(A/) 



il-f{M)fAUM) + 

dM.M.^{Ms,z)^M^ [z, Ms) 



(18) 



M 



'dMs 



IV. 



MAPPING THE GAMMA-RAY 
ANNIHILATION SIGNAL 

A. Galactic contributions 



The gamma-ray fiux from local substructures receives 
contributions from all the subhalos along the line of sight, 
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typically ^ 10^ subhalos when integrating over a a solid 
angle ^ 10~^sr. A brute-force integration with a Monte 
Carlo approach is therefore impossible even on high- 
speed computers. To circumvent this problem, Fieri et 
al. in Ref. [l^ have proposed a hybrid approach that 
consists in splitting the integral into two different con- 
tributions. The first one, which we regard as due to un- 
resolved subhalos, is the average contribution of a sub- 
halo population distributed according to Eq. [5l which 
can be estimated analytically. The second one comes 
from the nearest subhalos, that one may hope to resolve 
as individual structures, and it is estimated by numer- 
ical integration of 10 independent Monte Carlo realiza- 
tions. To determine the number of individual subhalos 
in each Monte Carlo realization we assume that their 
contribution to the gamma-ray flux represents a Poisson 
fluctuation to the mean gamma-ray annihilation signal. 
Following this criterion, we only consider, in each mass 
decade, subhalos with $'=°>5™° > ('I>'g'™° ^ (V' = 180°)) - 

10^^ GeV^ cm~^ kpc sr, where the brackets indicate the 
mean annihilation flux. For a given halo mass, this re- 
quirement corresponds to generate all subhalos within a 
maximum distance dmax{M). More than 500 subhalos 
are found within dmax for each mass decade from 10 AIq 
to IO^Mq and their positions are simulated according to 
Eq. O For less massive subhalos, if less than 500 subha- 
los are found, we increase dmax to include all the nearest 
500 subhalos in that mass range. 

Fieri et al. in Ref. [l^ have checked that the predicted 
gamma-ray flux does not depend on dmax ■ Here we need 
to perform a similar robustness test to check that the 
choice of dmax does not introduce spurious features in 
the angular power spectrum. For this purpose we have 
generated different Monte Carlo realizations simulating 
subhalos up to n x dmax with n from 2 to 5. We have 
checked that convergence in the power spectra is obtained 
already for n ~ 2. Since the convergence depends weakly 
on the model concentration parameter adopted, in this 
paper we have adopted a more stringent criterion and 
set our maximum simulation distance at 3 x dmaxiM) to 
guarantee that the convergence is reached in both sub- 
halo models considered. 

To summarize, in our model we consider three sources 
of gamma-ray annihilation within the MW: the smooth 
Galactic halo (which we will refer to as NFW), the un- 
resolved subhalos (UNRES), and the resolved substruc- 
tures (RES). Since the total annihilation depends on the 
square of the DM density, for the purpose of computing 
the angular correlation signal we also need to consider 
the double products NFW x RES and NFW x UNRES. 
We neglect RES x UNRES since a resolved subhalo at 
a given location excludes the presence of unresolved sub- 
structures. 

We point out that the subhalos labelled as "resolved" 
do not correspond to those substructures that will be 
detected by Fermi LAT. Instead they merely represent 
the nearest subhalos that we have generated to account 
for the discreteness of the subhalo spatial distribution. 



Moreover, considering the population of RES subhalos 
also accounts for the pixel-to-pixel variation in the an- 
nihilation flux that would be neglected focusing only on 
the smooth UNRES signal. 

The possibility of detecting galactic DM substructures 
is discussed in details in Ref. [l^. Since, in model 
Bza,ref, Fermi LAT is expected to detect only very few 
substructures, in this work all the simulated subhalos are 
included in the RES map. On the contrary, for model 
Bzf,ref, Fermi LAT may resolved up to ~ 100 point-like 
sources. Since we are interested in the correlation prop- 
erties of the diffuse gamma-ray emission, these sources 
should be masked. However, we checked that the angu- 
lar spectrum for the model B^f^ref does not significantly 
depend on the exclusion of these point-like sources. 

In Figures [T] and [2] we show the different contributions 
to the differential flux at 10 GeV as a function of the angle 
^ from the GC in model Bz„^ref and Bzf^ref, respectively. 

For Bza,ref (Fig. [1]) the flux from galactic substruc- 
tures (red curve) dominates over the smooth NFW pro- 
file (black curve) at angular distances from the GC larger 
than about 3° and it remains the dominant contribution 
at larger angles. The blue hne labeled 'EGRET' rep- 
resents the EGB inferred from EGRET measurements, 
equal to 1.10 X lO^^^cm-^g-iGeV^^sr-i at 10 GeV, that 
can be used as a generous upper bound (the Fermi LAT 
satellite will impose more stringent constraints, see next 
section). 

Note that the NFW curve in Fig. [T] exceeds the EGRET 
EGB in the central degree, but not the flux of the bright 
gamma-ray source (J1746-2851) that has been identified 
by EGRET to be near - possibly coincident with - the 
GC, and it is thus consistent with observations. 

For Bzj ,re/ , galactic substructures dominate the NFW 
signal at even smaller angular separations form the GC. 
In this case, the average flux for galactic substructures is 
larger than in the B^^^ref case and almost matches the 
EGRET background. 



B. Extra-galactic contribution 

The green line in Figures [1] and [2] represents the 
average extra-galactic flux contributed by DM annihi- 
lation in extra-galactic halos, and their substructures, 
calculated with Eq. [T31 In model i?zo,re/ this flux 
turns out to be 2.20 x 10^^^cm~^s~^GeV~^sr^^, signif- 
icantly smaller than in model Bz^.ref (for which it is 
1.45 X lO-i^cm-^s-iGeV^sr-i) but stiU wefl below the 
galactic contribution and the EGRET constraint. 

We stress the fact that in this work we have considered 
the EGRET background since it constitutes the best ob- 
servational constraint available to date. This background 
is likely to be dominated by unresolved extra-galactic 
sources like blazars. Using the blazars Gamma-ray Lu- 
minosity Function derived from the blazars detected by 
EGRET, Ando et al. have shown that these sources 
should contribute only to 25-50% of the total EGB [ll|. 



7 



NFW 

Galactic substructures 
Extragalactic halos and subhalos 
EGRET 




10' 



10-^ 



NFW 

Galactic substructures 
Extragalactic halos and subhalos 
EGRET 




J I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

20 40 60 80 100 120 140 160 180 
[degrees] 

FIG. 1: Different contributions to the differential flux 
d^/dEdO, [cm-^s"^GeV"^sr-i] at 10 GeV as a function of 
the angle 'I' from the GG, in model -BzQ.re/- The blue line cor- 
responds to the EGRET estimate of the EGB as parametrized 
in Ref. [l^. The contribution from DM annihilation in the 
smooth host halo and its substructures are represented by 
the black and red curves, respectively. The latter is ob- 
tained by averaging over 10 different Monte Carlo realiza- 
tions of the subhalos population. Bumps and wiggles are due 
to the contribution of the individual subhalos of the RES 
population. The green line represents the extra-galactic flux 
contributed by DM annihilation within extra-galactic halos 
and their substrutures. We have assumed = 40 GeV, 
av = 3 X 10~'^^cm'^s~^ and annihilations to bb. 

One therefore expects that, thanks to its superior sen- 
sitivity, Fermi LAT wiU be able to resolve a significant 
fraction of these sources [3^. As a result, the amplitude 
of the unresolved gamma-ray background will decrease, 
effectively lowering the blue line in Figs. [T] and [2] in the 
case of Fermi LAT. This will have three advantages. First 
of all, a fainter gamma-ray background will increase the 
probability of detecting the gamma-ray annihilation sig- 
nal produced by individual structures. Second, the back- 
ground from Fermi LAT will provide tighter constraints 
to our model for the DM annihilation, telling us if our de- 
scription is too optimistic. Finally, with a larger number 
of resolved astrophysical sources, Fermi LAT will provide 
a better estimate of the blazar luminosity function. 

The fact that the extra-galactic flux is contributed by 
both DM annihilation and blazars and because of the un- 
certainty on the latter contribution, one has the freedom 
of modifying the relative importance of blazars emission 
to the total EGB, keeping a good agreement with the 
data. Different choices for these relative contributions 
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FIG. 2: Different contributions to the differential flux 
d^/dEdO, at 10 GeV as a function of the angle from the GC 
in model Bzf,ref- The curves represents the same quantities 
shown in Fig. [T]but for -Bz/.re/- 



have been proposed in Ref. [131]. They are listed in the 
first two columns of Tab. [TTl The same DM annihilation 
flux leads to a different value of /nS""""^ for EGRET and 
Fermi LAT, due to the lower Fermi LAT EGB. 

It is particularly important to note that not all of the 
scenarios in Tab. [Til ^^re physically plausible when im- 
plemented in our model. For example, boosting the av- 
erage galactic and extra-galactic DM annihilation flux 
so that it will account for more than 61% of the Fermi 
LAT EGB, leads to a flux that exceeds the EGRET con- 
straint towards the GC (even excluding the central de- 
gree). Thus, for Bzg^ref we will only consider f^"^"" as 
large as 0.61, restricting to the last two rows in Tab. [Til 
For the model Bzf.ref, a similar argument rules out cases 
with ■ > 0.80. 



V. ANGULAR POWER SPECTRUM OF THE 
GAMMA-RAY UNRESOLVED SIGNAL 

A. Galactic contribution 

To compute the galactic gamma-ray angular power 
spectrum and investigate the relative importance of the 
different contributions, we separately analize the five 
maps of corresponding to the NFW, RES, UN- 

RES, NFW X RES and NFWxUNRES contributions. In 
the maps, the value of $cosmo specified within angular 
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0.1 


0.9 


0.03 


0.97 


ruled out 


ruled out 


0.3 


0.7 


0.20 


0.80 


ruled out 


5.9 


0.5 


0.5 


0.39 


0.61 


49.8 


4.2 


0.7 


0.3 


0.61 


0.39 


29.9 


2.6 



TABLE II: f^i'i^f^ and /oS"'"'^ indicate the contribution of 
gamma-ray from blazars and from DM annihilation respec- 
tively, to the EGB as measured by EGRET, /u™, and fS^m"' 
represent the same quantities estimated for the Fermi LAT 
satellite in Ref. [l^l using the blazar luminosity function de- 
rived from EGRET data. The last two rowes indicates the 
corresponding boost factors needed to bring the average an- 
nihation flux (galactic plus extra-galactic) to the relative per- 
centage of the EGB value. 



bins of An = 9.57 x lO^^sr. 

Since we arc interested in the contributions of DM ha- 
los and subhalos, we mask out the region close to the 
GC and the Galactic plane in which the signal is dom- 
inated by gamma-rays produced by cosmic rays inter- 
acting with the interstellar galactic medium [lj|. This 
background rapidly decreases with galactic latitude. For 
this reason we have used a composite mask consisting 
in a strip of 10° above and below the Galactic plane 
and the squared area around the GC with coordinates 
\b\ < 30° and |Z| < 30°. Outside the mask, the average 
galactic annihilation flux turns out to be (d^/dEdfl) — 

"^sr^^ in the case of Bza,ref 
1.16 X lO-^cm-^s-iGeV"^' 



and (d^/dEdil) ^ 1.16 x lO-^cm-^s-iGeV-'sr^i for 
Bzf,ref, and it is dominated by the UNRES term. 

To characterize the angular correlation signatures of 
the various contributions we compute the angular power 
spectrum of the different maps with the HEALPix 2.01 
package [13, EH- Five Healpix_Map objects are created 
covering the whole sky with 786432 pixels of constant 
area (corresponding to N_side = 2®). Each pixel of 
the Healpix_Map is filled with the corresponding flux 
d^/dEdO, from the simulated maps. The number of pix- 
els is determined from the requirement that the area of 
the single bin is smaller than Afi. 

The angular power spectrum, Cf, is computed from 
the spherical harmonic coefficients „i of the gamma- 
ray flux angular fluctuation map as follows: 



/ d$ \ 



\ dEdn I 



dn f'"^'^ 



/ d^ \ 



dEdn \ dEdn / 



2^+1 



(19) 



(20) 



where Yi^m{0^ (j)) are the spherical harmonics and the in- 
tegration in Eq. [12] extends to the unmasked area. The 
angular spectra of the RES and NFWxRES maps have 
been obtained by averaging the spectra of the 10 inde- 
pendent maps corresponding to the different Monte Carlo 
realizations of resolved subhalos. 



We note in passing that the angular correlation signal 
does not depend on the particle physics scenario since 
it is evaluated with respect to the average flux (aside 
from a weak dependence on the minimum subhalo mass, 
which in turn depends on the properties of the DM parti- 
cle). However, for the purpose of comparing the gamma- 
ray angular spectrum due to DM annihilation to that of 
extra-galactic astrophysical sources, a particular particle 
physcis scenario has to be specified. In this case we refer 
to the reference values presented in Section IIII Al 

Our masking procedure induces spurious features in 
the NEW, UNRES and NFWx UNRES spectra that we 
smoothed out by averaging over 18 multipoles. 

The total angular spectrum is obtained by performing 
a sum of the aforementioned contributions and the cross- 
correlation terms: 



dEdn / ^ 



T.{dEdn) 



(21) 



^\ dEdn / \ dEdn / 



with i,j e [NEW, RES, UNRES, NEWxRES, 
NFWx UNRES]. 

The results for model Bzg^ref are plotted in Fig. [31 
for which we also have taken into account the pixel 
window function for the corresponding resolution of our 
maps [4]| . The different curves in the plot repre- 
sent the angular spectrum of each component weighted 
by the square of its relative contribution to the total 
annihilation flux {d^' / dEdn)'^ / {d^^°' / dEdnf . Cross- 
correlations are taken into account in the computation 
of the total spectrum but are not shown in Figs. [3| and 
[H Black crosses represent the spectrum of the total flux 
and the shaded areas show the Icr error boxes. These are 
obtained by summing in quadrature the uncertainties re- 
lated to the finite binsize and, if applicable, the scatter 
among the different Monte Carlo realizations of resolved 
subhalos. The first contribution dominates at lower mul- 
tipoles I while the second becomes the main source of 
error at larger multipoles where the power spectrum is 
mainly due to the RES contribution. Poisson errors are 
not accounted for since they depend on the observational 
setup and they will be discussed in the next Section for 
the case of the Fermi LAT satellite. 

The main contribution to the total flux is provided by 
'the autocorrelation of the resolved (red curve) and unre- 
solved (blue curve) subhalos. The former dominates at 
large multipoles while the latter dominates at f < 30. 
The interplay between these two components is respon- 
sible for the minimum at £ ~ 30 in the total spectrum. 
All components involving the smooth DM distribution 
provide a negligible contribution to the total spectrum. 

Fig. [H shows the angular spectrum of model B^j^ref- 
The total spectrum and the different contributions are 
very similar to model Bzg^ref- The only difference is rep- 
resented by the position of the minimum which is now 
around i ~ 50. This is due to the contribution of the 
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FIG. 3: The curves represent the angular power spectrum 
Ci£{£+ 1) /2n of the different contributions to the galactic an- 
nihilation flux in the framework of model B^g^ref- The labels 
identify the different contributions. Black crosses represent 
the angular spectrum of the total annihilation signal together 
with its la error (shaded areas). 

resolved subhalos (red curve) which is smaller than in 
the previous case: subhalos in the B^f^ref niodel are less 
concentrated than in Bzo.ref and consequently they con- 
tribute less to the spectrum at small angular scales. We 
conclude that changing the concentration parameter sig- 
nificantly affects the mean annihilation flux but not the 
angular power spectrum of the signal. 

Our model predictions depend on the assumption 
about the minimum subhalo mass, for which the the- 
ory docs not provide very strong constraints. In order 
to check the dependence of our results on the minimum 
subhalo mass, we have re-simulated the resolved subha- 
los excluding all structures below a given mass threshold. 
The spectrum does not change when excluding subhalos 
with a mass below W^Mq, showing that the angular cor- 
relation is mainly contributed by massive subhalos. 

To check the reliability of our hybrid model, we have 
compared our results with those of a similar study re- 
cently carried out by Siegal-Gaskins [HI who used a pure 
Monte Carlo approach to compute the power spectrum 
of the annihilation signal produced by a population of 
subhalos in the MW above lOA/© and IO'^Mq, hence ne- 
glecting the contribution of smaller subhalos. To com- 
pare our simulation approach to hers, we have modified 
our subhalo model to match that of Ref. [2l|. In par- 
ticular we have re-simulated our galactic subhalo using 
a shallower mass function dn/dM oc M~^'^, assuming 
an Einasto density profile phaios for subhalos, adopted a 



FIG. 4: The curves represent the angular power spectrum 
Cil{l -f l)/27r of the different contributions to the galactic 
annihilation flux in the framework of model -B^^. ,re/. 



simplified power-law model for the concentration param- 
eter with c(M) (X M~^-^^^ and normalized the number 
of subhalos in the MW with a mass larger that lO^M© 
to be N{M > IQ^Mq) = O.OO64(1O8M0/Mmw)"''-^ (this 
normalization forces the number of MW subhalos to be 
five times smaller than with the normalization assumed 
in Sec. [n]for the B^g^ref and B^^^ref cases). 

The only residual difference between our model and 
that of Fig. 5 and 6 of Ref. [2l| is the spatial subhalo 
distribution in the MW which is more similar to the anti- 
biased model of Ref. [2l| . 

The result of this exercise is shown in Fig. [51 The 
curves shows the spectra normalized at = 150 for the 
various cases explored. The spectra of the subhalos a la 
Sieg al-Gaskins are consistent with Fig. 5 and 6 of Ref. 
[2l| for both the two mass thresholds considered (10 AIq 
(black curve) and 10''Mq (red curve)). The agreement is 
valid at the same time for the shape and the amplitude 
of the signals. 

The comparison with the spectrum of our subhalos in 
the case of Bza^refi shows that, for th e s ame mass cut, 
our spectrum is steeper than in Ref. [2l|. This can be 
explained by noticin g th at small mass halos are more 
concentrated in Ref. [2l| than in our case, resulting in a 
larger power at small scales (i.e. at large multipoles). 

Also the amplitude of the angular spectra predicted in 
models B^g^ref and B^^^ref results to be different than 
the same quantity in Fig. 5 of Ref. [2l|, even when 
the different subhalo spatial distribution is considered: 
our spectra are characterized by a lower amplitude, the 
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FIG. 5: Angular power spectrum Ct, normalized to Ciso for 
the resolved subhalos modeled as in Ref. [2ll]. The black 
(red) curve shows the case of a IOMq (10'^ Mq) mass cutoff. 
The green (blue) curve indicates the spectrum for the resolved 
subhalos in model B^^^ref for a mass cutoff of IOMq (10^ Af©). 



result of having considered subhalos with masses lower 
than 10 A/0. They contribute to the total annihilation 
flux but their distribution is almost homogeneous around 
the observer, hence decreasing the spectral amplitude. 

Finally, as we have already pointed out, the results 
presented here are robust to the choice of the maximal 
simulation distance. We have checked that resolving ha- 
los at distances larger than 3 x dmax{M) has very little 
impact on the angular power spectrum. 



B. Extra-galactic contribution 

Our analysis for the angular spectrum of the extra- 
galactic annihilation flux is similar to the one originally 
developed by Ando et al. in Ref. [l3| and extended 
in to account for substructures. Like in our model, 
these authors assume that DM halos and subhalos have 
NFW density profiles and that subhalos trace the smooth 
underlying mass distribution of the host halo. Moreover 
they assume that the mean number of subhalos within a 
host halo of mass M, {N\M), scales as oc and they 
consider two extreme cases, a = 1 and a = 0.7 [isj . In 
this work we only consider the first case to match the 
subhalo mass function of our model. 

We notice that Ando et al. [l^l do not account for the 
contribution of subhalos by integrating over their mass 
function, as in the second term of Eq. [181 Instead, they 
consider a gamma-ray subhalo spectrum averaged over 



its mass function and multiplied by the mean comoving 
density of subhalos. These model details hardly affect 
the angular power spectrum since one considers fiuctu- 
ations about the mean flux. While we will quantify the 
differences between the two approaches in a future work 
in which the angular power spectrum and the mean flux 
will be computed self-consistently, here, for the purpose 
of discussing the angular power spectrum, we adopt the 
same approach as in Ref. and we will use their very 
same angular power spectra shown in Fig. 5 of their 
paper for both the DM annihilation and the blazars con- 
tribution to the extra-galactic flux. 

We do not plot the contribution of the DM and blazars 
to the angular spectrum of the extra-galactic gamma-ray 
background in Figs. [3]and|4]to avoid confusion. Instead, 
in the next Section we show how well these spectra will 
be measured by Fermi LAT and compare them to the 
galactic contributions. 



VI. FERMI LAT GAMMA-RAY ANGULAR 
SPECTRA 

In this section we apply our model to predict the an- 
gular power spectrum that Fermi LAT is expected to ob- 
serve and discuss the possibility of indirect DM detection 
through characteristic signatures in the spectrum. 

The contribution of the different components to the 
total angular spectrum of the gamma-ray flux is shown 
in Fig. [6] for the i3zo,re/ case. The red histogram indi- 
cates the angular spectrum of anisotropics in the blazar 
emissivity (see Fig. 5 of Ref. [l^). Black crosses rep- 
resents the contribution of galactic substructures that, 
as we have shown, dominate over the smooth NFW sig- 
nal. Open diamonds indicate the contribution of extra- 
galactic halos and their substructures. All the different 
contributions are normalized to the square of the ratio 
between the average flux of that particular contribution 
to the average total flux {d<i>ydEdfl)'^/{d^""/dEdnf. 

Poissonian errors in each bin of the sky maps are due to 
the small statistics, i.e. the limited number of gamma-ray 
photons from annihilation collected during observations. 
The la error is computed, following Ref. [ij], as follows: 



5C, 



{21 + i)a;/,„.. 



Ci + cf + 



c 



N 



(22) 



where A/ = 18 and 47r/f „^ = 9.706 sr is the area out- 
side the mask. Cjv is the power spectrum of the photon 
noise Cat = 47r/(„^ /TVegb that depends on the instru- 
ment characteristics and integration time (A'^egb is the 
number of photons of the EGB) and is independent from 
We have assumed an effective area constant with en- 
ergy of lO^cm^ and an exposure time of 1 year. Wi is 
the window function of a Gaussian point spread function 
Wi — exp(— ^^(7^/2) that we compute for the Fermi LAT 
angular resolution ~ 0.115°. The angular spectrum of 
the background Cp is not uniquely defined but depends 
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on the signal one wants to detect: for the extra-galactic 
DM component, the background is represented by blazars 
and by the galactic annihilations. While, viccversa, for 
the galactic DM component, blazars and extra-galactic 
annihilations play the role of background. 

The total signal (shown as crosses in Fig. [5] 
and Fig. [7]) is obtained by adding all components 
weighted by their relative contribution to the total 
gamma-ray flux {d<^^ /dEdn)'^ / {d<P^°ydEdn)^ including 
the cross-correlation term involving extra-galactic halos 
and blazars. We have ignored the cross-correlation be- 
tween the galactic subhalos and the blazars since they 
have independent spatial distributions. The total errors 
(blue boxes) account for the Poisson noise of both the 
galactic and extra-galactic DM component, to which we 
summed the binsize and scatter among Monte Carlo real- 
izations. The Poissonian noise becomes the main source 
of errors only al large multipoles. Finally, we have as- 
sumed that the angular power spectrum of blazars is 
known without errors. 

In the particle physics scenario considered in this paper 
and for the Bzo,ref model, the average flux {d^ / dEdH.) 
produced by DM annihilations within and outside our 
galaxy is less than 1% of the EGB estimated for Fermi 
LAT. We can therefore boost up this contribution by in- 
creasing the cross section of the DM candidate. This has 
the effect of changing the values of f^i^^l^ and /o^^J"', i.e. 
the relative contribution of blazars and DM to the total 
EGB. The two panels in Figure [6] refer to f^Z"" = 0-61 
(left) and /d™" = 0-39 (right). The corresponding boost- 
ing factors are shown in the plots. Larger boosting fac- 
tors (corresponding to larger /om'"0 are excluded since 
they would correspond to models already excluded by 
the current EGRET constraint. 

In both plots the contribution of the galactic signal 
to the angular power spectrum largely dominates the 
extra-galactic component at all multipoles. As a conse- 
quence, the two main contributions to the angular spec- 
trum are provided by the blazars and the galactic annihi- 
lation signal. The former dominates at large multipoles 
whereas the latter dominates at low i. The position of 
the crossover depends on the boosting factor, but even 
in the less favorable case of /^m"' = 0-39 the DM annihi- 
lation signature can be clearly seen as a turnover in the 
power spectrum at ^ < 30. 

These considerations remain valid for the Bzf.ref 
model whose angular power spectrum is shown in Fig. 
[71 In this case, however, the crossover is found at smaller 
multipoles, making more difficult to detect the contri- 
bution for DM annihilation when /^m"' = 0.39 because 
of the large errorbars. On the other hand, with the 
Bzf.ref model one can increase the boost factor up to 
/dm"" = 0-80 without ending up with an unphysical sce- 
nario already excluded by EGRET. 



VII. DISCUSSION AND CONCLUSIONS 

In this paper we have modeled for the first time the 
angular power spectrum of the diffuse gamma-ray signal 
at high galactic latitude taking into account both the DM 
annihilation in substructures within our Galaxy and an- 
nihilation within extra-galactic halos and their subhalos, 
also including the contribution from extra-galactic ob- 
jects like blazars. The purpose is to understand whether 
the gamma-ray satellite Fermi LAT will be able to unam- 
biguously detect the signature of DM annihilation that, 
in the ACDM scenario, preferentially takes place in the 
central regions of DM halos and in their hierarchically- 
nested substructures. 

For this purpose we have developed an hybrid approach 
in which the annihilation signal is computed using two 
different techniques. First, wc perform a numerical inte- 
gration to compute the contributions to the gamma-ray 
flux from the smooth galactic halo, from the unresolved 
galactic substructures and, finally, from extra-galactic 
DM halos and their substructures. Second, we use Monte 
Carlo techniques to account for the nearest galactic sub- 
halos, that show up in the mock gamma-ray maps as 
individual gamma-ray sources. 

We have explored two different subhalo models corre- 
sponding to different prescriptions for the mass concen- 
tration within the halo, Bz„.ref and B^f^ref- 

The EGB is computed by adding the annihilation sig- 
nal from extra-galactic halos and that from high energy 
astrophysical sources like blazars. The relative contribu- 
tions of the two signals is treated as a free parameter. 
The only constraint is that the total gamma-ray signal, 
galactic annihilation foreground and EGB, does not ex- 
ceed the current limit of EGRET. 

Our main conclusions are: 

• The annihilation of the smooth galactic halo dom- 
inates over that of galactic substructures only 
within a few degrees from the GC. The exact value 
depends on the subhalos concentration and ranges 
from < 1° in model Bz^^ref, to ^ 3° in model 
Bzo,ref- Moreover, the gamma-ray flux from extra- 
galactic halos is fainter than the galactic signal at 
all angles. This is true both for Bzg.ref and Bzj^ref- 

• To predict the annihilation flux we have assumed a 
favourable particle physics setup, with a mass of 40 
GeV, a cross section a?; = 3 x 10~^^ cm~'^s~^, and 
annihilations to hb. With this choice, model Bzo.ref 
accounts for only 0.8% of the total unresolved EGB 
that will be measured by Fermi LAT. Larger fluxes 
are possible for smaller masses and larger cross sec- 
tions, provided that the total flux does not exceed 
the current measurement of the EGB by EGRET 
and still allowing blazars to contribute to a signif- 
icant fraction of the flux. Similar considerations 
also apply to model Bz^^ref- The largest boost fac- 
tor compatible with the data is of order 100, and it 
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FIG. 6: Angular power spectrum of gamma-ray anisotropies in the case that DM annihilations contribute to 61% (left) and 
39%(right) of the total average flux of the EGB as estimated for the Fermi LAT satellite. The plus signs + indicate the 
galactic component following the BzQ.re/ model (see Fig. [3]). The diamonds o refer to DM annihilating in extra-galactic, 
unresolved DM halos and subhalos and the spectrum is taken from Ref. [l^. The red line represents the angular spectrum 
of blazars as described in Ref. [ij. All these three components are normalized to their average contribution to the total flux 
(d$ydEdQ,f/{d$*^°^/dEdnf. The Icr error of the total spectrum is indicated by the blue boxes and it is obtained propagating 
the errors from the galactic and the extragalactic DM components, assuming no errors for the blazars. 
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FIG. 7: Angular power spectrum of gamma-ray anisotropies in the case that DM annihilations contribute to 80% (left), 
61%(center) and 39% (right) of the total average flux of the EGB as estimated for the Fermi LAT satellite. The curves and 
the error boxes have the same meaning of Fig. |6]but refer here to the B^^^ref case. 



could be achieved e.g. through the so-called Som- 
merfcld enhancement or in some rather fine-tuned 
supersymmetric scenarios [13, 113 • 

• The angular power spectrum predicted by our 
model is reliable and robust. To test its reliability 
we compared our results with those recently pub- 
lished in Ref. [2l| , who performed a similar analysis 
considering, however, only galactic substructures 
and using a full Monte Carlo approach. When the 



same halo model is considered, we find fair agree- 
ment with Ref. (211 . To test the robustness of our 



method we have verified that our results do not 
change significantly when we increase the volume 
within which we simulate the distribution of galac- 
tic subhalos beyond a well defined reference value. 

• The angular power spectrum of the galactic 
gamma-ray photons produced by DM annihilation 
is dominated by subhalos. Contributions involv- 
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ing the DM particles in the smooth galactic halo 
are negligible at all multipoles. At low multipoles 
(large angles) the spectrum is dominated by un- 
resolved halos whose annihilation flux appear in 
the gamma-ray sky map as a smooth foreground 
(see e.g. Figs. 6 and 7 of Ref. 0). On the 
other hand, nearby resolved clumps dominate the 
power at small angular separations (large multi- 
poles). These results are robust in the sense that 
they do not change appreciably when adopting dif- 
ferent prescriptions for the concentration parame- 
ters. This is not surprising since the power spec- 
trum, which refers to fluctuations about the mean 
flux, should not depend appreciably on the subhalo 
structure. 

Subhalos of masses below IO^Mq give a small con- 
tribution to the angular spectrum of the galactic 
annihilation flux. 

The signature of DM annihilations in the angu- 
lar power spectrum can be found only for an op- 
timistic particle physics setup. In fact, for our 
benchmark particle physics scenario (without any 
boost factor) the power spectrum of the unresolved 
gamma-ray background would be completely domi- 
nated by blazars. However, boosting up the particle 
physics factor without exceeding the background 
that EGRET has estimated, we find that the power 
spectrum is dominated by the DM component at 
low multipoles while blazars determine the spec- 
trum at small angular scales. The turn-over de- 
pends on the prescription for the concentration of 
subhalos, occuring at smaller scales in the case of 

Bzo,ref- 

The angular power spectrum of DM annihilations is 
largely dominated, at all scales, by the galactic sig- 
nal. It depends mostly on the relative importance 
of the average galactic flux compared to the extra- 
galactic DM signal which, in our models, favours 
the first term, both for Bz^ ^ef and Bzj.ref- 

We find that a 1-year all-sky survey with Fermi 
LAT, now operational, may be able to spot the an- 
nihilation signature of DM in the angular power 
spectrum of the unresolved gamma-ray flux. This 
signature is provided by the up-turn in the power 
spectrum which, at low multipoles, is dominated by 
the DM annihilation signal. Our results partially 
contradict those of Ref. [ll] who did not account 
for the galactic contribution to the angular power 
spectrum which, instead, completely obliterate the 
one provided by annihilations in extra-galactic ha- 
los. 

Our results seem to confirm those of Ref. [2lj 
since the angular spectrum measured by Fermi LAT 
should indeed provide an indirect detection from 
DM annihilation within our Galaxy, although it 



is doubtful that such measurement will be able to 
constrain the properties of the subhalo population. 
However, we should point out that the subhalo de- 
scription adopted in Ref. [2l| would boost up the 
extra-galactic DM annihilation flux to a value of 
3.60 X 10~^"cm"^s^^GeV^^ comparable to, if not 
larger than, the galactic annihilation flux. On one 
side this fact increases the possibility of detecting 
the DM contribution to the angular spectrum. On 
the other hand it will be hard to disentangle galac- 
tic and extra-galactic contributions. 

As a remark we stress that the results of the recent Via 
Lactea II [lOl and the Aquarius [3l|, [12 numerical experi- 
ments, both with improved resolution, have now updated 
our knowledge on the subhalo distribution function and 
concentration parameters. The differences in the subha- 
los extracted from these simulations, whose main char- 
acteristics arc listed in Table HI may change the angular 
spectrum of the DM annihilation flux although we have 
shown that changing subhalo model significantly affects 
the gamma-ray flux but has less impact on the spectral 
analysis. We will investigate this issue thoroughly in a 
forthcoming paper. Here we just report the result of a 
preliminary analysis in which we have computed the an- 
gular power spectrum of the galactic annihilation signal 
obtained when subhalos are modeled according to Via 
Lactea II and Aquarius. In both cases we flnd that the 
angular spectrum is flatter, especially at small I. This is 
not surprising since, in both experiments, the mass distri- 
bution is less clumpy and the subhalo radial distribution 
is flatter than in the model adopted in the present paper. 
A flatter spectrum, combined with a fainter annihilation 
flux due to the reduced number of substructures, would 
reduce the up-turn in the angular spectrum, making it 
more challenging to detect the DM annihilation feature. 
The signiflcancc of this effect depends on the boosting 
factors that, because of the dimmer flux, could be fur- 
ther increased without exceeding the measured gamma- 
ray background. 

Finally, being tuned to the results of N-body simula- 
tions, our description of the galactic and extra-galactic 
subhalo population docs not account for the effect of 
baryons in the evolution of DM structures. As a conse- 
quence our computation of the angular spectrum does not 
account for the effect, e.g., of Intermediate Mass Black 
Holes, which may be present in our Galaxy and which 
correlation properties have been recently studied in Ref. 
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